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Abstract: Superfluid condensates are known to occur in contexts rang- 
ing fi'om laboratory liquid helium to neutron stars, and are also likely to 
occur in cosmological phenomena such as axion fields. In the zero tem- 
perature limit, such condensates are describable at a mesoscopic level by 
irrotational configurations of simple relativistic perfect fluid models. The 
general mechanical properties of such models are presented here in an in- 
troductory review giving special attention to the dynamics of vorticity flux 
2-surfaces and the action principles governing both individual flow traje- 
tories and the evolution of the system as a whole. Macroscopic rotation 
of such a condensate requires the presence of a lattice of quantised vortex 
defects, whose averaged tension violates perfect fluid isotropy. It is shown 
that for any equation of state (relating the mass density p to the pressure 
P) the mesoscopic perfect fluid model can be extended in a uniquely simple 
and natural manner to a corresponding macroscopic model (in a confor- 
mally covariant category) that represents the effects of the vortex flbration 
anisotropy. The limiting case of an individual vortex defect is shown to be 
describable by a ( "global" ) string type model with a variable tension T (ob- 
tained as a function of the background fluid density) whose "vorton" (i.e. 
closed loop equilibrium) states are derived as an exercise. 
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1 Introduction. 



These lectures offer an introduction to the dynamics of vortices in simple 
superfluid models of a general category that includes the kind appropri- 
ate in the context of laboratory condensed matter for the representation of 
Helium-4 in the zero temperature limit (for which no "normal" entropy car- 
rying constituent is present), as well as the ultrarelativistic "stiff" kind that 
has been widely used ^ for the representation of a massless axion field in 
cosmology. Between these extremes, this category also includes the kind of 
model suitable as a simple approximation for the description of the neutron- 
pair condensation superfluid that is generally believed (for both theoretical 
and observational reasons 0) to occur in the intermediate layers of neu- 
tron stars. (A more sophisticated treatment would include allowance for 
the presence of an interpetrating ionic crust lattice or of an independently 
superconducting proton pair condensate, not to mention complications such 
as spin.) 

The declared purpose of this school for which these lectures are intended 
is to study analogies between low temperature laboratory phenomena and 
high energy cosmological phenomena. However it is desirable, when possible, 
to reinforce mere analogy by interpolation. I would therefore like to empha- 
size the particular interest for this purpose of the intermediate regime of 
neutron star interiors, about which a great deal of admittedly indirect infor- 
mation is available from analysis of pulsar frequency variations. In contrast 
with most other areas of astrophysics, but in common with typical labora- 
tory applications of condensed matter physics, the temperatures in typical 
neutron star interiors can usually be considered to be very low compared 
with the relevant energy scales. On the other hand, in contrast with typical 
laboratory applications of condensed matter physics, but in common with 
high energy cosmological scenarios, the treatment of neutron star interiors 
requires allowance for significant relativistic effects. 

While evidently indispensible for an accurate treatment of neutron star 
matter, a relativistic treatment is commonly - but wrongly - considered to 
be an unnecessary complication for low temperature laboratory applications. 
For a reasonably accurate treatment of laboratory liquid Helium a relativis- 
tic treatment is indeed unnecessary, but what is wrong is to suppose that 
it is more complicated. On the contrary, as I hope these lectures will make 
clear, the relativistic treatment, in so far as it is available, is mathematically 
simpler. (The essential reason for this is that the Lorentz group is, in the 
technical sense "semi-simple" whereas the Galilei group is not.) 
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The work, in collaboration with David Langlois j^, on which these lec- 
tures are based, was originally inspired by the observation by Rick Davis 
and Paul Shellard Q (in precisely the spirit that is the official raison d'etre 
for the present school) of a strong analogy between the behaviour of vor- 
tices in the "stiff" (ultrarelativistic) massless axion model at one extreme, 
and in the incompressible superfluid model that is commonly used for the 
description of zero temperature Helium-4 at the opposite extreme. 

What is shown here is the way to carry out the analogous, but not quite 
so trivially simple, treatment that is needed for the intermediate category 
of generically compressible superfluid models. At a mesoscopic level these 
models will be of irrotational perfect fluid type, and will be characterised by 
a subluminal speed, Cj say, of ordinary "first" sound, that will be determined 
~ by an equation of state specifying the pressure P as a function of the mass 
density p — according to the familiar formula 

c, = {dP/dpf'\ (1) 

Rotation at a macroscopic level will entail the presence of an Abrikosov 
type lattice of quantised vortex defects whose averaged tension produces a 
deviation from perfect fluid isotropy. Analysis of individual vortex tubes ^] 
indicates that, in the limit for which their relative "drift" or "flight" velocity 
is highly subsonic, their averaged effect can be described by the inclusion of 
an extra isotropy violating term of a uniquely simple and natural kind in 
the action for a corresponding macroscopic model. Despite their violation 
of isotropy, such models (for different equations of state) will be shown to 
belong to a category that is conserved by conformal transformations. More 
particularly, it will be shown that the specific model corresponding to the 
high pressure limit of a gas of negligibly interacting particles as characterised 
by c'^ = (where c is the speed of light) has the same kind of conformal 
invariance property as the well known example of Maxwell's equations. 

The scope of the present review does not extend to the recently per- 
fected relativistic analogue |^, |^ of Landau's original non-dissipative two 
constituent superfluid model (which has recently been shown to be express- 
ible in a very elegant Galilean covariant form |8|) for the purpose of allowing 
for the dynamical effect of the "normal" entropy flux current that would be 
present at a non zero temperature). Such effects are not the ones that are 
most important in the context of neutron stars, the main application for 
which the work described here is intended, where the temperature will typi- 
cally be so low (compared with the relevant M.e.v. range energy scales) that 
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thermal corrections can for most purposes be neglected. For application to 
neutron stars, a more important kind of generalisation is to a relativistic 
models capable of describing neutron superfluid penetration |^] of the solid 
material forming the crust, and of describing [^0| the protonic supercon- 
ductivity that is expected Q within the neutron superfluid at a deeper level 
below the crust. 



2 Canonical treatment of relativistic flow trajec- 
tories. 

Before proceeding it is desirable to recall some essentials the relativistic 
kinematics and dynamics, particularly in view of the regretable tradition in 
non-relativistic fluid theory - and most notably in non-relativistic superfluid 
theory - of obscuring the essential distinction between velocity (which for- 
mally belongs in a tangent bundle) and momentum (which formally belongs 
in a cotangent bundle) despite the fact that the distinction is generally re- 
spected in other branches of non-relativistic condensed matter theory, such 
as solid state physics, where the possibility of non-alignment between the 
3- velocity u", and the effective 3- momentum pa of an electron travelling in 
a metallic lattice is well known. 

In a non relativistic treatment it is only in strictly Cartesian (rather 
than e.g. cylindrical or comoving) that the distinction between contravari- 
ant entities such as the velocity v"" and covariant entities such as the momen- 
tum Pa can be ignored. In a relativistic treatment, even using coordinates 

<-> {t,x"} of Minkowski type, with a flat spacetime metric Qf^i, whose 
components are of the fixed standard form diag{— c^, 1, 1, 1}, (where c is 
the speed of light) the necessity of distinguishing between raised and lowered 
indices is inescapable. Thus for a trajectory parametrised by proper time r, 
the correponding unit tangent vector 

is automatically, by construction a contravariant vector: its space compo- 
nents, u"" = 7f" with 7 = (1 — u^/c'^)"^/^ will be unaffected by the index low- 
ering operation ^ u^j, = g^j.uu'^ but its time component = dt/dr = 7 
will differ in sign from the corresponding component uq = —^(? of the as- 
sociated covector u^. On the other hand the 3-momentum pa and energy 
E determine a 4- momentum covector vr^y that is intrinsically covariant, with 
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components Tr^ = Pa, t^o ~ ~^ ■ '^^^ covariant nature of the momentum 
can be seen from the way it is introduced by the defining equation, 

in terms of the relevant position and velocity dependent Lagrangian function 
L, from which the corresponding equation of motion is obtained in the well 
known form 

^ = ^ (4) 

In the case of a free particle trajectory, and more generally for fluid flow 
trajectories in all the simple "barotropic" perfect fluid models with which 
the present lectures will be concerned, the Lagrangian function will have the 
familiar standard form 

L = ^fJ-QfiiyU^u" - \^c^ , (5) 
in which (unlike what is needed for more complicated chemically inhomoge- 



neous models [11, ^]) it is the same scalar spacetime field /i that plays the 
role of mass in the first term and that provides the potential energy contri- 
bution in the second term. The momentum will thus be given by the simple 
proportionality relation 

7r„ = fiUjy , (6) 

so that one obtains the expressions E = Pa = fJ-JVa, in which the field 

H is interpretable as the relevant effective mass. 

In the case of a free particle model, the effective mass /x will of course 
just be a constant, ^ = m. This means that if, as we have been supposing 
so far, the metric g^y is that of flat Minkowski type, the resulting free 
particle trajectories will be obtainable trivially as straight lines. However 
the covariant form of the equations (|2|) to (^) means that they will still be 
valid for less trivial cases for which, instead of being fiat, the metric g^y 
is postulated to have a variable form in order to represent the effect of a 
gravitational field, such as that of a Kerr black hole (for which, as I showed 
in detail in a much earlier Les Houches school [^], the resulting non trivial 
geodesic equations still turn out to be exactly integrable). 

In the case of the simple perfect fluid models with which we shall be 
concerned here, the effective mass field ^ will be generically non-uniform. 
In these models the equation of state giving the pressure P as a function 
of the mass density p can most conveniently be specified by first giving p 
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in terms of the corresponding conserved number density n by an expression 
that will be decomposible in the form 

p = mn + 4t , (7) 

in which m is a fixed "rest mass" characterising the kind of particle (e.g. a 
Cooper type neutron pair) under consideration, while e represents an extra 
compression energy contribution. The pressure will then be obtainable using 
the well known formula 

P = (n/i - py , (8) 

in which the effective dynamical mass p (or equivalently the "specific en- 
thalpy" pc^) is given by 

dp _ ^ 1 de 
dn dn 

It is this that is to be taken as the effective mass function appearing in the 
specification (|5|) of the relevant Lagrangian. 

When one is dealing not just with a single particle trajectory but a 
spacefilling fluid flow, it is possible and for many purposes desirable to con- 
vert the Lagrangian dynamical equation (0) from particle evolution equa- 



= 3r = "^ + 72 — • (9) 



tion to equivalent field evolution equations [11, H]. Since the momentum 



covector vr^ will be obtained as a field over spacetime, it will have a well 
defined gradient tensor VpiTu that can be used to rewrite the left hand side 
of (^) in the form diTy/dT = uP\IpiTy. Since the value of the Lagrangian 
will also be obtained as a scalar spacetime field L, it will also have a well 
defined gradient which will evidently be given by an expression of the form 
VyL = dL/dx^ + {dL/d-Kp)VyT:p . We can thereby rewrite the Lagrangian 
dynamical equation (Q) as a field equation of the form 

UPVpTT^ + TTpV^UP = V^L . (10) 

An alternative approach is of course to start from the corresponding 
Hamiltonian function, as obtained in terms of the position and momentum 
variables (so that formally it should be considered as a function on the 
spacetime cotangent bundle) via the Legendre transformation 

H = TT^U''-L. (11) 

In this approach the velocity vector (^) and the dynamical equation (Q) are 
recovered using the familiar formulae 

(h^_dH_ dH_ 
dr diTu ' dr dx'^ ' 
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The consideration that we are concerned not just with a single trajectory 
but with a spacefilhng fluid means that, as in the case of the preceding La- 
grangian equations, so in a similar way this familiar Hamiltonian dynamical 
equation can also be converted to a field equation which takes the form 

2ufV[p7r,] = -V,H , (13) 

with the usual convention that square brackets are used to indicate index 
antisymmetrisation. On contraction with u'^ , the left hand side will evidently 
go out, leaving the condition 

u''V^H = 0, (14) 

expressing the conservation of the value of the Hamiltonian along the flow 
lines. 

The actual form of the Hamiltonian function that is obtained from the 
particularly simple kind of Lagrangian function (|5|) with which we are con- 
cerned will evidently be given by 

H = —g-^PTT^TTp + ^ . (15) 

In order to ensure the proper time normalisation for the parameter r the 
equations of motion (in whichever of the four equivalent forms (^, ([lOD, 
(|l2|) , ([l3[) may be preferred) are to be solved subject to the constraint that 
- in order for to be correctly normalised - the numerical value of the 
Hamiltonian should vanish, 

H = =^ uf'Uf, = -c^ , (16) 

initially , and hence also by (|l^) at all other times. In the more general 
systems that are needed for some purposes the Hamiltonian may be con- 
strained in a non uniform manner [|ll|, ^ so that the term on the right of 
(|l|) will be non zero, but in the simpler systems that suffice for our present 
purpose the restraint (|l6|) ensures that this final term will drop out, leaving 
a Hamiltonian equation of the very elegant and convenient form 

^^^V[^7r,] =0. (17) 
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3 Vorticity conservation and flux 2-surfaces. 



The preceding form (^) of the dynamical equations is particularly handy for 
the analysis of symmetries and the derivation of conservation laws according 
to which physically interesting quantities are preserved by various kinds of 
continuous displacement Q . The variation induced by an infinitesimal 
displacement generated by an arbitrary vector field say, will be given by 
the corresponding Lie derivative operator, whose effect on a scalar field, n 
say, will evidently be given simply by 

kCn = k''Vufi, (18) 

while its effect on the metric will be given by the well known though not 
quite so immediately obvious formula 

kCg,,, = 2V(^/c,) , (19) 

using round brackets to denote index antisymmetrisation. If k^ is a Killing 
vector field, i.e. if the displacement generated by k^ is a symmetry of the 
spacetime metric, then the right hand side of (jl^) will vanish. In curved 
space such symmetries are rare, but in ordinary fiat space timethere is of 
course a ten parameter family of such Killing vector fields generating the 
Poincare group, whose algebra has as its basis the four independent gener- 
ators of uniform spacetime translations and the six independent generators 
of the Lorentz group. The effect of the Lie differentiation operation on an- 
other vector field, u^, say will be given simply by their mutual commutator 
bracket, 

kCuP = [k, u]P = -uCkP , [k, uf = k^V^u^ - k^V^k^ . (20) 

Such Lie differentiation is an example of an operation for which the dis- 
tinction between covectors and ordinary contravariant vectors is important: 
except when performed with respect to a Killing vector. Lie differentiation 
does not commute with index raising and lowering. The rule that applies to 
the covector = g^^u^ can evidently be obtained by combining (|^) and 
( pO| ) , so that for -n^ one obtains a formula that can be conveniently expressed 
in terms of exterior (antisymmetrised) derivatives in the form 

kCiT^ = 2FV[p7r^] + V^(7rpF) . (21) 

Although it does not commute with index raising and lowering. Lie differen- 
tiation does commute with exterior differentiation. Thus for the relativistic 
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vorticity tensor, which is defined to be the antisymmetrised derivative of the 
momentum covector, i.e. 

Wfj,i, = 2V^7r^] , (22) 
so that its own exterior derivative will automatically vanish, i.e. 

V[^u;,p] = 0, (23) 

it follows that its Lie derivative will be obtainable just by taking the exterior 
derivative of (|2^) which gives 

kCw^^ = -2Vi^{w^]pkP) . (24) 

The preceding general formula (^) can immediately be used to rewrite 



the Lagrangian dynamical equation (|10[) in the expressive form |11, 



uCtt^ = VyL , (25) 

from which by taking the exterior derivative, one can immediately derive the 
dynamical vorticity conservation law that is the key to superfluidity theory 
(and much else) in the form 

uCw^y = . (26) 

The interpretation of this crucially important result is that for any fiow 
governed by the Lagrangian equations ^ or their Hamiltonian equivalent 
(|T2|) (which have so far simply been postulated ex cathededra, but whose 
validity for any "barotropic" perfect fiuid be made clear in the following 
section) the vorticity field ( p2|) will simply be convected onto itself by the fiow 
field m'^, with the implication that if it vanishes initially w^y will remain zero 
throughout the flow, which in this case will be describable as "irrotational" . 

The foregoing results can be considerably strengthenned in cases such as 
those of the simple "barotropic" perfect fluids considered here for which the 
proper time normalisation is ensured by the Hamiltonian constraint ([T^. 
This has the effect of reducing the dynamical equations to the particularly 
simple form ([T7|) , whose interpretation is that the fiow vector u^^ must be 
an zero eigenvalue eigenvector of the vorticity tensor w^^- The posession of 
a zero eigenvalue requires that Wf^y should satisfy the degeneracy condition 

w^[uWp„i^ = , (27) 
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which excludes the possibility of it having matrix rank 4, with the impli- 
cation that unless it actually vanishes it must have rank 2 (since an anti- 
symmetric tensor can never have odd integer rank). This means that the 
flow vector is just a particular case within a whole 2-dimensional tangent 
subspace of eigenvectors satisfying 

e^Wf,^ = . (28) 

This subspace will be spanned by a unit worldsheet element tangent bivec- 



tor £f^f^ of the kind whose use was developped by Stachel |13], and that is 
definable, wherever the vorticity magnitude 

w = C^w^.w^^n'^'' (29) 
does not vanish, as being proportional to the dual vorticity tensor W^'^, i.e. 

w 2 



gt^'^ = - W^\ W^" = -e'^^"" Wp^, (30) 



(note that sign convention used for the worldsheet orientation here is the 
opposite of what was used in the preceding article (P) where ef^'^^"' is the 
totally antisymmetric tensor normalised by the convention that its non zero 
components are equal to 1 or -1 (depending on whether the index ordering 
is an even or odd permutation) with respect to locally Minkowskian coor- 
dinates with gfj,u diag{— c^, 0, 0, 0}, so that with respect to an arbitrary 
coordinate system the non zero component values of gi^'^P^ will be given by 
±c\\g\\-^/\ 

It can be seen that this worldsheet tangent bi-vector will satisfy 

£^"'£,,u = -2c\ £^'•'w,p = 0. (31) 

It follows from this last equation that the contraction of any covector with 
£^^'^ will provide a solution of the vortex worldsheet tangentiality condition 
(p8|). A noteworthy example is the helicity /i'^ as defined |[ll|, ^ by 

hf" = w£^''7r^ , (32) 



which can be seen from (^) and ( p7D to satisfy the helicity current conser- 
vation law [0, 

V^h" = . (33) 

Again using the degeneracy condition (^) and the condition that the 
vorticity also satisfies the Poincare closure condition (p3|) , it can be shown M 
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that the tangent elements characterised by £^'^ and generated by solutions 
of (^) (or more specifically by vy and hy) will automatically satisfy the 
relevant Frobenius condition for integrability, meaning that they will mesh 
together to form well behaved timelike 2-dimensional worldsheets. This 
makes it possible to extend the flow line conservation laws resulting from 
any continuous symmetries that may be present. 

The simplest example is the Bernouilli type theorem that applies - even 
if the Hamiltonian does not satisfy the constraint (^) - whenever is a 
symmetry generator of the system, so that in particular kCH and kCiTy 
both vanish: it can be seen from from and that u''Vu{T^pkP) wih 
vanish, so TTpk^ will be constant along each flow line. In the most obvious 
application, the Killing vector is just the generator of time translations in 
flat space with the Minkowski coordinates, k^ <-> {1, 0, 0, 0}, so that the 
Bernouilli constant will simply be identifiable with the negative of the effec- 
tive energy per particle, i.e we shall have 'KpkP = —E with E = which 
automatically includes allowance for both compression energy and kinetic 
energy contributions. 

This conclusion can be greatly strengthenned when the Hamiltonian con- 
straint (|16|) holds: one then gets e'^'V^i'n'pkf) = for any vector e'^ satis- 
fying (|2^) which evidently means that TTpk'^ will not just be constant along 
each flow trajectory but that it will be constant throughout each of the 
2-dimensional vorticity flux worldsheets. In the irrotational case 

wp^ = 0, (34) 

for which (^) is satisfied trivially by any vector e'^ at all, one can draw the 
even stronger conclusion that iTpk^ will be constant throughout the fluid: 

kC7rp = V^(7rpF)=0. (35) 

(A similar conclusion of global uniformity of TTpk^ would also be obtainable 
immediately from ( p!6|) and (|2l| ) in the alternative, more widely familiar, case 
of motion that is rigid in the sense of having a flow vector that is aligned 
with the Killing vector, i.e. for which u^'^k'^^ = 0.) It is to be emphasised 
that the applicability of this kind of generalised Bernouilli theorem is not 
limited to the case of ordinary stationarity, for which k^ is a time translation 
generator, but is just as well applicable to cases of axisymmetry for which 
the Killing vector is a rotation generator, and it has recently been found 
to be very useful |jl5| for a hybrid case (involving a non-rigidly rotating 
binary pair of tidally deformed and thus non-axisymmetric neutron stars 
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which, from the point of view of a distant observer, are not stationary but 
periodically evolving, but which are nevertheless stationary from the point 
of view of a local observer with respect to a suitably rotating frame). 

An alternative way of obtaining the local vorticity conservation theorem 
(p6|) is as the differential limit of the global Kelvin-Helmoltz theorem to the 



effect that, as is manifest from the form (25) of the Lagrangian dynamical 
equation, the action integral 

■Kydx^ (36) 

will be preserved if taken round a closed circuit that is convected by the 
flow field u^. In the irrotational case (^), the Jacobi type action integral 
(|36|) between any two fixed endpoints will be unaffected by continuous dis- 
placements of the path between them, and hence can be used to construct 
a locally well defined field S that is not only such that one has 

T^u = VuS , (37) 

but that will also be a solution of the Hamilton Jacobi equation specified by 
setting the Hamiltonian function to zero with the gradient of S substituted 
in place of the momentum. 

The special case of a (simple, zero temperature limit) superfiuid is speci- 
fied by the existence of a well defined mesoscopic phase factor e^^ (represent- 
ing the phase factor of an underlying bosonic condensate that might consist 
of Helium-4 atoms or Cooper type neutron pairs) in which the phase angle 
If is given according to the usual correspondence principle by 

ip = S/n. (38) 

In a multiconnected configuration of a classical irrotational fluid the Jacobi 
action fleld S obtained from ( |37| ) might have an arbitrary periodicity, but 
in a superfiuid there will be a U(l) quantisation requirement that the peri- 
odicity of the phase angle (p should be a multiple of 2tt, and thus that the 
periodicity of the Jacobi action S should be a multiple of 2Trh. The simplest 
configuration for any such superfiuid is a uniform stationary state in a flat 
Minkowski background, for which the phase will have the standard plane 
wave form 

S/h = kax" - ojt , (39) 

from which one obtains the correspondence vTjy <-> {—fioj, hka}, which means 
that the effective energy per particle will be given by £" = Ji^c^ = fiu) and 
that the 3- momentum will be given by pa = n^Va = hka- 
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It is to be remarked that for ordinary timelike superfluid particle tra- 
jectories the corresponding phase speed oo/k oi the wave characterised by 
( |39| ) will always be superluminal, a fact which people working with liquid 
Helium-4 in the laboratory can blithely ignore, since what matters for most 
practical purposes is not the phase speed but the group velocity of perturba- 
tion wave packets. Our present discussion will be limited to the strict zero 
temperature limit for which no such packets are excited, but it is easy to 
extend the relativistic analysis to low but non-zero temperatures for which 
the relevant excitations are phonons 0. Although their phase speed and 
group velocity are the same, both being given by the formula (jl]) for the 
ordinary ("first") soundspeed which will of course be subluminal, phonons 
do nevertheless have a tachyonic aspect of their own: their 4-momentum 
covector is always spacelike, in contrast with that of an ordinary fluid or 
superfluid particle which is timelike. This means that whereas the effective 
energy E of an ordinary fluid or superfluid particle is always positive, the 
effective energy E oi a, phonon may be positive or negative, depending on 
whether the frame of reference with respect to which it is measured is mov- 
ing subsonically or supersonically. The well known implication is that if the 
superfluid is in contact with a supersonically moving boundary there will 
inevitably be an instability giving rise to dissipative phonon creation. 

4 Conventional formulation of perfect fluid and 
simple superfluid theory. 

Although sufficient for the derivation of many important properties of the 
flow, the dynamical equations on which the preceding sections have been 
based contain only part of the information needed for a complete determi- 
nation of the perfect fluid evolution. The most usual way of presenting the 
complete set of equations of motion of a simple perfect fluid of the barotropic 
type we are considering - meaning one whose intrinsic local physical state 
is characterised just by a single independent scalar field - is in the form of 
a conservation law of the standard form 

V^T^'' = , (40) 

for a stress momentum energy density tensor that is specified as a function 
of the timelike unit flow tangent vector and a single independent scalar 
field variable, such as the conserved particle number density n, on which the 
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other relevant quantities, such as the effective mass /i given by will be 
functionally dependent. 

For a system of this simple kind, the 4 independent components of the 
energy momentum conservation law ( ^0|) provide all that is needed to deter- 
mine the evolution of the 4 independent components that characterise the 
local state of the system, which can be taken to be the scalar n and the 3 
space components of (since the remaining component u° is not dy- 
namically independent but determined by the unit normalisation condition 
( [l6|) as an algebraic function of the 3 other components). 

In any perfect fluid model the mass density p and pressure P are physi- 
cally characterised by their role in the specification of the stress momentum 
energy density tensor, for which the standard expression is 

T^"' = (p+ ^^u^'u" + Pg^"' . (41) 

In the simple barotropic case, the relation (P) between the dependences of 
P and p on n can be seen to be necessitated by the requirement that the 
dynamical system (^) should ensure conservation of the number current 

= nu'' . (42) 

It is easy to check, using (^) that contraction of (^) with does indeed 
lead to the required result, namely 

V^n'^ = . (43) 

One can also verify the not quite so well known result |14| that the remaining 



independent equations (|40| ) can be reorganised in the canonical uniformly 
Hamiltonian form ( [T7| ) on which the work |jll|, ^ of the preceding subsection 
was based, and which is expressible succinctly as 

n^Wf,^ = . (44) 

Thus in addition to the formula (|4^ ) that has been used so far, the only 
additional information needed for the complete specification of the dynamics 
of a barotropic fluid system is the obvious particle conservation law (|4^). 

Given a dynamical system, one of the first things any physicist is inclined 
to ask is whether it is derivable from a Lagrangian type variation principle. 
We have already seen in the previous sections that (|4^ by itself is obtainable 
from Lagrangian equations of motion for the individual trajectories, which 
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are of course obtainable from a one dimensional action integral of the form 
J Ldr with L as given by ^ . The question to be adressed now is how to 
obtain the complete set of dynamical equations ( [40D , including (^) as well 
as (^4|), from an action integral over the 4-dimensional background manifold 
5'*' of the form 

J = y £ dS^''' , dS^'^ = ^^^^d'^x , (45) 

for some suitable scalar Lagrangian functional C. 

Several radically different procedures are available for doing this. Al- 
though ultimately equivalent "on shell", they involve variation over "off 
shell" bundles that differ not just in structure but even in dimension. The 
oldest and most economical from a dimensional point of view is the worldline 



variation procedure developed by Taubfl^], followed Clebsch type variation 
procedure developed by Schutzp^, but for our present purpose it will be 
more convenient to employ the more recently developed Kalb-Ramond type 



method! 18| that has been specifically designed for dealing with problems of 
macroscopic superfluidity. 

The problem is greatly simplified if, to start off with, one restricts one- 
self to the purely irrotational case (^), which is all that is needed for the 
description of zero temperature superfluidity at a mesoscopic level. For this 
case independent variable can be taken to be just the Jacobi action S, or 
equivalently in a superfluid context, the phase (p as given by (|^), and the 
action is simply taken to be the pressure P expressed as a function of the ef- 
fective mass n, with the latter constructed as proportional to the amplitude 
of the 4-momentum, according to the prescription 

//2c2 = (46) 



with the 4-momentum itself given by the relation ( p7\ ) that applies in the 
irrotational case, i.e. 

vr^ = nVuif . (47) 

Thus setting 

C = P, (48) 

and using the standard pressure variation formula 5P = c^nSfi one sees that 
the required variation of the Lagrangian will be given by 

6C = -rfdiTy = -hn''V^{6ip) . (49) 
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Demanding that the action integral ( |45| ) be invariant with respect to in- 
finitesimal variations of then evidently leads to the required conservation 
law (|]). 



5 Introduction of the dilatonic amplitude field 



For an equation of state such that P oc /i^, the Lagrangian ( ^8|) will, as it 
stands, have the quadratic field gradient dependence that is typical of simple 
physical field theories. This occurs in the special limit case, relevant to the 
massless axion field in cosmology fl], ^, of the "stiff" model, characterised 
by ( |62[ ) as discussed below, for which the "first" sound speed (|T|) is equal to 
that of light. Except in this idealised limit feature of the preceding 

Lagrangian function ( p8| ) that is widely considered to be undesirable is its 
generically non quadratic dependence on Vip, an apparent drawback that is 
not uncommonly dealt with by recourse to approximation [|l9| ]. 

What this secion will show however is that that, not just for P (x fj? but 
for a quite general equation of state, the Lagrangian (^) can be reformulated 
in the much more desirably fashionable form 

^ = -f ^'(V.¥^)V>-yW (50) 

- in which the potential energy density term V is some suitably chosen 
algebraic function of the amplitude $ - by a transformation of variables 
that is absolutely exact [18|, without any need for recourse to approximation 



provided one adopts the correct definition for the auxiliary field variable 

What would involve an approximation would be to to implement a fur- 
ther step whereby $ is "promoted" from the status of an auxiliary variable 
to that of an extra dynamical variable by adding in a supplementary kinetic 
term, 

A£{,| = -a27i2(v-$)V,$, (51) 

for some constant value of a. Various kinds of gradient term, of which this 
is the most obvious, might indeed be added for the purpose of improving 
the physical precision of the model in regimes of rapid amplitude variation. 



where deviations from the strictly isotropic perfect fiuid form (41) might 
be expected to become significant. However for a strongly interacting liq- 
uid like Helium-4 (as opposed to a weakly interacting gas) I know of no 
reliable procedure for the quantitative evaluation of such a term on the ba- 
sis of an underlying many-particle quantum theory. The most commonly 
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used [19 1 ansatz, namely to take a = 1, is not automatically guaranteed 
to provide an improvement but might even bring about a deterioration of 
physical precision in some circumstances. The standard choice a = 1 does 
not have theoretical or empirical foundations but is merely based on the 
purely mathematical consideration that it provides an adjusted Lagrangian 
Ci^iy = C + ^{1}^ of the "relativistic Ginzburg Landau" form 

Ai} = -f (52) 

in which the phase variable ip and the amplitude variable <I> have been com- 
bined to form a complex variable 

$ = 6^"^, ^ = e-''^, $ = 1^-1. (53) 

This feature is useful for some purposes ~ notably for providing a smoothed 



out treatment [pO| , 21] of the vortex defects that arise where the phase (p 
becomes indeterminant - but not so convenient for deriving exact results 
such as the uniformity of the Bernouilli constant ir^k'^ that was shown above 
to be valid for a model of the original kind whenever k'^ is a symmetry 
generator. 

As far as physical accuracy (as opposed to mathematical convenience) is 
concerned it will usually not matter very much whether one uses the exact 
perfect fluid model given by ( |50| ) or the associated Ginzburg Landau model 
given by (52) since their difference, as given by the extra kinetic term (51) 



can be expected to be very small except in the immediate neighbourhood of 
a vortex core where neither model can be expected to be physically accurate. 
What does matter for physical accuracy is the choice of the functional form 
for the amplitude 

When it comes to the point, without even attempting to provide any 
serious microscopic derivation, many introductory presentations start by 
assuming the validitity of the Landau Ginzburg framework, and then resort 
to crude guesswork for the specification of <I>, typically taking <I> tx y/n, or 
even more commonly 

^>oc^. (54) 

Such a choice just happens to provide a quantitatively acceptable result in 
the non-relativistic limit, but the flimsiness of its theoretical basis is exposed 
by the fact that for a generic relativistic model it gives the wrong answer 
in the regime of slow amplitude variation where the (correctly calibrated) 



Landau Ginzburg model (52) should agree with the relevant perfect fluid 
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model ( |50|). The correct choice - the only way to get exact agreement 
between ([48D and (|5^) - is to take 



^ /n\i/2 



the corresponding characterisation for the potential energy density function 
V being [^] that it should be given by 

V='-^. (50) 



When characterised in this manner, - with the correct identification (|55| ) 



rather than (54) - the perfect fluid model given by the Lagrangian (^) pro- 
vides a uniquely canonical model for the representation at zero temperature 
of a simple superfluid such as Helium-4, at least in the regime where the 
rate of variation of the amplitude $ is not too rapid compared with that 
of the phase. It is to be emphasised that for this regime the justification 
for this perfect fluid model is actually on a sounder footing than that of the 
corresponding Ginzburg Landau type model ([5^) , since appart from the pos- 



tulate of the isotropic perfect fluid form (41) for the stress momentum enery 
density tensor, whose applicability in the slow amplitude variation regime is 
hard to doubt, and the invocation of the generally valid conservation laws 



(HO) and (42), the only assumption on which this fluid model is based, and 



for which it relies on an underlying microscopic quantum analysis, is that of 
the rather well established existence of the mesoscopic phase ^p. 

Having evaluated y as a function of ^> one can recover the effective mass 
/i, number density n, mass density p and pressure P of the fluid using the 
formulae 

^^=Za^^^ n = $V, (57) 



and 



p=\ + ^ , P = \ ^V^c^ - V , (58) 

which are derivable from (^) and (|^). 

Since the preceding relations entail the variation rule 

5V = i?c^^5<^, (59) 

it can be seen that the full variation of the reformulated perfect fluid La- 
grangian ( |50|) will be given by 

5C = -{ttuTt" + p^c^)<^6^ - h^^TT^V^idip) , (60) 
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with the momentum covector tTi, as specified before by (p7|). Thus in this 
new formulation, instead of imposing (^) as a defining relation, we obtain 
it as a field equation from the requirement that the action integral should 
be invariant with respect to local variations of the auxiliary field while 
as before the analogous requirement for variations of the phase variable (p 
gives back the particle conservation law (^2|) . 



5.1 Special equations of state 

It is to be remarked that whereas for a generic compressible fluid equation 
of state giving p as a function of n, and hence giving P as a function of 
/i, the formula ( |59|) will provide a corresponding function V that provides 
a reformulated action function of the form ( |50| ) from which the required 
perfect fluid dynamical equations are obtainable by treating $ and if as 
independent. There is however an exceptional case that works somewhat 
differently, namely that of the "stiff" Zel'dovich model - pertaining to the 
massless axion field in cosmology ^ - which is characterised by 

P = pc\ (61) 

corresponding by (Q) to sound propagation at the speed of light. This case is 
obtained from a primary equation of state of the form p ccv? oi equivalently 
(as observed at the beginning of Section |) P oc /i , which gives y = 
and p oc n. This last relation means that the amplitude ^ will simply be a 
constant, and so will not act as an auxiliary variable in the usual way. In this 
case what one has to do to obtain the relevant "stiff" dynamical equations 
is simply to treat if as the only independent variable in the Lagrangian (|50| ) . 
There is another important special limit for which it is not the new 



formulation but the original variational formulation (48) in terms of the 
pressure function that fails to work, namely the extreme case - which often 
useful as an approximation - of a "dust type" model characterised by p oc n 
for, which there is no pressure, i.e. P = 0. This low pressure limit model 



(the only model for which the naive identification (54) is exactly valid) is 
immediately obtainable within the framework of the new formulation (|50| ) 
simply by taking V (x 

The simplest non trivial potential energy function that can be used in 
the new formulation ( |50[ ) is provided by the "radiation gas" model (|6l| ) 
characterised by 



c2 



c, 



' = V ^ yoc$^ (62) 
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which is obtainable from a primary equation of state of the form p oc n"^/^. 
This is the model that is appropriate for representing the cosmological black 
body radiation (which is to a good approximation irrotational, though of 
course it is not a superfluid). It also applies to a degenerate Fermi gas of 
massless (or due to high compression effectively massless) non interacting 
particles, and for that reason was used as a first crude approximation in some 
of the pioneering studies of neutrons stars, whose intermediate layers are 
indeed believed to be superfluid. As will be explained below, this particular 
model is characterised by conformal invariance of the same kind as is familiar 
in the well known case of Maxwellian electromagnetism. 

For more general purposes, as an approximation that can be usually be 
expected to be reasonably accurate within a limited density range, one can 
of course use the obvious generalisation of (^) to the standard form 

for suitably adjusted constants m and a. For such a model the fluid mass 
density and pressure will be given according to (^8|) by 



p = ^2^2^3^^4^ P = a2$\ (64) 



from which it can be seen that the sound speed (the quantity one would 
probably want to use in practice to fix the appropriate value of a) will be 
given according to (||) by 



2 /3 \-i 



< = (::2 + ^] (65) 



6 Introduction of the Kalb Ramond gauge field 

The treatment given in the preceding section is ideal for the description of a 
superfluid at a mesoscopic (i.e. intervortex) scale, but for the treatment of 
an ordinary perfect fluid with rotation, or for the treatment of a superfluid 
on a macroscopic scale (allowing for the averaged effect of a large number 
of vortices) more general models are required. 

A first step towards the kind of generalisation that is needed is to for- 
mulate the current in terms of an antisymmetric Kalb Ramond type tensor 
field B^i, = —By^ whose exterior derivative 

N^,p = 3V[^i?,^] , (66) 
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ie^-^'^iV,,. (67) 



whose (Hodge type) dual 



is to be identified with the particle number current (p^D, which will evidently 
be invariant under the effect of Kalb Ramond gauge transformations Bfj_i, 
Bf^i, + 2V[nXu] ■ It can then be seen that the consequent closure condition, 

V[^iV.p.]=0, (68) 



is equivalent to the usual form (|43| ) of the particle conservation law. 

The idea now is to perform a Legendre type transformation C ^ A 
whereby the independent scalar field ip of the preceding formulation based 
on C is replaced by the antisymmetric gauge tensor i?^,^ in a new formulation 
of the same model in terms of a different dually related Lagrangian function 
A which takes the form 

^ = -r^N^'''N,,p-Vm. (69) 



In the preceding formulation based on C as given by ( pOD the irrotationality 
property was kinematically imposed in advance while the particle current 
conservation was obtained from the variational principle as a dynamical 
equation. However in the reformulated version based on the dual Lagrangian 
( |69| ) it is the particle current conservation law that is obtained in advance 
as we have seen via the kinematic identity (|68[), while on the other hand the 
irrotationality condition (|3^), i.e. 



V[,.^.]=0, (70) 

is obtained directly from ( p9| ) in the equivalent dual form 

V^i'^-'^N'^P'') =0. (71) 

from the requirement of invariance with respect to independent variations 
of the gauge 2-form B^^. 

The purpose of replacing the simple scalar field (p by the tensorial field 
Bfj^i, is to enable the extension of the model to the general perfect fiuid case, 
in which the particle conservation law (^) is retained, but the irrotationality 
condition (^) is abandoned. The way to do this |jl^ is to introduce a 
complete Lagrangian of the form 

£ = A - ie^"''" 5^.10^. , (72) 
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where, in analogy with what has already been done for the current 3-form, 
the vorticity 2-form is constructed from independent gauge fields in such a 
way that its conservation property, (^3|) is automatically ensured in advance 
as a kinematic identity. For this purpose the requisite independent gauge 
fields can be taken to be a pair of independent scalars, say, in terms of 
which an identically conserved vorticity flux will be given by 

= 2(V[,x+)V,]X" ■ (73) 

When the Lagrangian (^) is substituted in the action integral (p5D, 
the requirement of invariance with respect to local variations of the (non- 
physical) dynamical gauge fields B^^, and x^, and of the (physical) auxiliary 
amplitude can be seen to lead back to our original dynamical momentum 
transport equation (^^. Thus (since the particle conservation law ( ^3[ ) has 
been imposed kinematically in advance) it provides the complete system of 
generic perfect fluid equations of motion as given by the standard stress 
momentum energy conservation law (^). 



7 Macroscopic allowance for vortex quantisation. 

All we have done so far is to reformulate ordinary barotropic perfect fluid 
theory in such a way that the vorticity w^i, comes in as an independent 
dynamical variable determined by the pair of scalar gauge fields x^- The 
model describing simple zero temperature superfluidity is obtained within 
this framework simply by taking the scalars to be constants (e.g. zero) 
so as to get Wij,^, = 0. 

For a macroscopic treatment of a bulk superfluid fibred by a dense con- 
gruence of discrete vortex defects, the need for an extra term in the action 
to allow for the extra tension and energy in the vortices was recognised long 



ago [22 1 in the laboratory context of Helium-4, and has more recently been 
taken into account in the Newtonian mechanical analysis of neutron star 
matter In an analogous manner, following earlier work |24] on the 



corresponding relativistic formulation needed for a more accurate treatment 
of neutron star matter, the more complete treatment ^ to be described 
here takes account of the averaged effect of quantised vortices aligned in an 
Abrikosov type lattice by summing over contributions of individual vortex 
cells as estimated Q using the usual (very good) approximation in which the 
hexagonal cells are treated as if they were cylindrically symmetric. What 
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this analysis suggests is that the averaged effect of such vortices can be rep- 
resented rather well - provided any relative flow is highly subsonic - by a 
remarkably simple and mathematically elegant modification of the generic 
potential function introduced by (|5^): all that seems to be necessary is to 
make an adjustment of the very simple form 

V{^} ^ + <^^T{w} , (74) 

for some function T depending just on the scalar magnitude w of the vor- 
ticity as given by (|2^). Furthermore, as was observed in analogous investi- 
gations in a non relativistic framework [^] , the dependence on the vorticity 
magnitude is approximately linear, having the form 

= Xw , (75) 

in which 

A = /C$2 , T = }Cw, (76) 

with a coefficient fC that can be taken to be a constant of the order of the 
Planck value, /C ~ ?i. 

A remarkable consequence of the ansatz ( |7^ in conjunction with the 
linearity postulate (76) is that, as will be described below, the extended 



model retains the noteworthy, though little known, conformal convariance 
property of the simple perfect fluid model. 

A more specific quantitative estimate of the value of the coefficient in 
([76|) is given by 

/C = ^, ^=ln{^}, (77) 

where 6 is an inner cut off length representing the microscopic vortex core 
radius, and A is a long range cut of length typically representing the mean 
intervortex separation distance. 

For a very precise treatment one would of course need to allow for a weak 
logarithmic dependence of /C on since will be inversely proportional 
to w, but so long as A is very large compared with 6, as will be the case 
in typical macroscopic applications, the effect of such a refinement will in 
practice be negligible. To obtain very high precision when the Mach (flow 
to sound speed) ratio is non-negligible [|2^ one might also have to allow for 
some sort of tensorial, not just scalar dependence on the vorticity tensor. 

Assuming that a sufficiently accurate treatment is obtainable without 
the need to take account of tensorial vorticity dependence, it follows that 
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the complete macroscopic superfluid Lagrangian will be given by 

^ = -4^2^'"'^^^^ - + - B,u)W^'' , (78) 



where W^'^ is the dual vorticity vector as defined by (pOl). 

In order to obtain the equations of motion for this system, one needs to 
evaluate the corresponding variation which will be expressible by 

6L = TT.dn-' - ^w^.db^" - '^{b^" + X^n^w^u + - ^ - (79) 

using the notation 

= ie^^/^^Sp^ , X"" = -w^"" . (80) 
^ w 

Requiring invariance of the corresponding integral with respect to SB^i,, or 
equivalently to 6b'^'^ , just leads back again to the usual relation (]2^) speci- 
fying the vorticity 2-form w^j^i, as the exterior derivative of the momentum 
1-form TTiy, so there will be a conserved helicity vector given by exactly the 
same formula (^) as before. The equation of state relation (^^ is how- 
ever modified by the addition of an extra term proportional to the vorticity: 
requiring invariance with respect to 6^ gives 

$V = -i:^ = ^-TT + 2Aw; . 81 

The final requirement is invariance with respect to variations of the gauge 
scalars determining the vorticity field according to (^), which is equiva- 
lent to the gauge independent requirement of invariance with respect to Lie 



transportation, Swf_ii, = kCw^y, as given by (p4D for an arbitrary displace- 
ment vector field k^^. The resulting dynamical equation is expressible in the 
form 

^^w^^ = , (82) 

which differs from its analogue in the perfect fluid limit only by the replace- 
ment of the conserved particle current by an "augmented" current 'k'^ 
that is also automatically conserved 

= , (83) 

whose specification is given by 

ttriA' = n'^ + VA^"" , (84) 
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which is equivalent, by ( pOD , to the condition that it be obtainable from a 
corresponding "augmented" Kalb-Ramond gauge 2-form given by 



(85) 



The geometrical interpretation of the modified dynamical equation (^) is 
facilitated by the observation that it is the same as the single vortex equation 



of motion (128) that will be derived below. 



8 The conformal covariance property 

It is well known that - when expressed in the usual way in terms of an 
electromagnetic gauge potential and charge current = en^, where e is a 
charge coupling constant and is an automatically conserved current vector 
such as is obtainable from a corresponding closed 3-form N^^p according to 
the Hodge type duality formula (^) - the equations of ordinary Maxwellian 
electromagnetism are preserved by any conformal transformation of the form 

9tMu ^ 9iiiy = e^'^9iiu , (86) 

for an arbitrary scalar field (j), on the understanding that the transformation 
affects neither the gauge 1-form, which obeys Aiy i— > Aiy, nor the closed 
current 3-form, which obeys N^^p i-^ Np^p. This last condition means that, 
since ( ^6| ) implies 

^pi^pa ^ -pupa ^ ^-i^^pvpa ^ ^g^) 

the current vector itself will undergo a conformal transformation of the form 



In the framework of the Kalb-Ramond representation (^6|) this is evidently 
equivalent to the requirement of preservation of the gauge 2-form, 

Bpv ^ Bpi, = Bpi, . (89) 

What is not so well known (since the utility of the auxiliary field as 
correctly defined by (^), is not yet widely appreciated) is that the ordinary 
barotropic perfect fluid equations (^) and (|4^) are also preserved by such a 
conformal transformation, subject to the understanding that the auxiliary 
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amplitude field undergoes a corresponding conformal transformation of the 
form 

q>^^ = e-^^ , (90) 
which is what is required to ensure the preservation, 

TTu ^ T^u = T^u , (91) 

of the momentum 1-form vr^ (whose role in the perfect fluid case is analogous 
to that of the gauge 1-form in the electromagnetic case). 

The preservation of the form of the dynamical equations (|4^ and (44) 



is in general not quite sufficient for preservation of the complete system, be- 
cause it is also necessary to satisfy the algebraical equation of state relation 
( |57[ ) specified by the function which governs the relation between the 

amplitude ^ and the 4-momentum magnitude fi = {—■k'^tTpY/'^. However 
provided the conformal scalar (j) is chosen to depend only the amplitude 

making the later a function of new variable <1>, then the system will be 
formally covariant in the sense that the new system will also behave as an or- 
dinary barotropic fluid but with a modified equation of state function 
in place of the orginal potental 

My original discussion of the perfect fluid case [^] envisaged a sce- 
nario involving gravitational coupling in the framework of a Brans-Dicke- 
Jordan type generalisation of Einstein's theory, with the dimensionless field 
(j) = ln{ <!>/$} acting as a dilatonic coupling scalar. Conformal covariance of 
this scheme was found to require that the transformation law for V should 
be given by an algebraic relation of the form V = V/^^'^. (For readers 
interested in this scenario I should warn that it is necessary to the correct 
the final sentence, concerning a special case for which there is a transforma- 
tion to a form in which "the dilatonic field is genuinely absent": the term 
"genuinely" should be replaced by "apparently", since although it disap- 
pears from the fiuid sector of the Lagrangian it effectively turns up again in 
the gravitational sector instead.) 

In present discussion I wish to describe a kind of conformal covariance 
like that of the well known Maxwellian example, having nothing to do with 
any particular kind of gravitational coupling theory, whether it be that of 
Einstein or anyone else. In order to preserve the formal structure of the fluid 
system by itself (without involving anything to do with active gravitational 
coupling) it can be seen from (|57|) that since we shall have 



^H^/i = e-V, (92) 
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it is necessary and sufficient that the effective potential should transform in 
such a way as to satisfy the differential condition 



^ — - = • (93 

d^ d^ ' 

This means, for example, that if the new system is to be characterised by 
a pressure free "dust" type equation of state of the form V = im^c^$^ for 
some constant rh (which has the convenient feature that the flow trajectories 
will simply be geodesies with respect to the new metric g^i,) then the desired 
transformation giving $ as a function of will be obtainable immediately 
by taking rh^cP^^^ = ^^"^dV/d^ except in the exceptional case for which 
the right hand side of ( p3| ) is constant. 

Whereas a perfect fluid solution for a generic equation of state can thus 
be conformally transformed to a solution for any other generic equation of 
state, there is an exceptional case for which such a procedure fails, namely 
the case with a constant value for right hand side of (^), which is that of 
the "radiation gas" model for which the dependence on $ of the effective 
potential V is of the homogeneously quartic form (|6^). It is apparent that 
the particular form 

2 

P = ^ oc $ (94) 

of the equation of state relation for this particular model will automatically 
be preserved by ( |90|) and (|9^), not just when the conformal scalar is chosen 
as a function of $ but for any field (p whatsoever. 

In summary, just as solutions of Maxwell's equation are well known to be 
mapped conformally onto other solutions of Maxwell's equations, solutions 
of the equations of the "radiation gas" model are analogously mapped con- 
formally onto other solutions of the "radiation gas" model, while generically 
solutions of the barotropic perfect fluid equations for any equation of state 
can be conformally mapped onto solutions for another different equation of 
state. In view of the intrinsic isotropy of a perfect fluid these conformal 
properties are not very startling. However we conclude this subsection by 
the rather more surprising observation that these properties will still hold, 
subject to the same condition (|9^), when the intrinsic isotropy is violated 
by the inclusion of the term (^) that allows for the energy and tension in 
quantised vortices. Since we shall have 
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and 

w w = e~'^'^w , (96) 

it can be seen that in order to have invariance of ^B^,y as given by ( |85|) and 
of the extra term Xw in the action integral (remembering that the measure 
dS^'^^ will transform proportionally to e^"^) we must have 

A^A = e~2'^A. (97) 

Rather remarkably, this condition just happens to be satisfied by the formula 
([76|), but it would fail for a more complicated vorticity dependent term such 
as might be needed for an accurate treatment of relative flow at a speed 
comparable with that of sound. 



9 The thin vortex string Hmit 

Let us now consider the limiting case for which instead of being 4- 
dimensionally extended, the vorticity distribution is concentrated in the 
neighbourhood of some particular vorticity flux 2-surface 5^^^ say, which 
might conveniently be characterised by zero values for the scalar gauge fields 
X^- Such a 2-surface will be describable in terms of internal coordinates o"" 
and fj^, which might conveniently be specified by the values on the world- 
sheet of the spacetime coordinates t and . 

However the internal coordinates may be chosen, the worldsheet em- 
bedding {cr", a^} = x^{a} will induce (in the technical sense "pull 
back" ) a two surface metric with components 'jab on the worldsheet that will 
be given by 

using a comma to denote partial differentiation with respect to the internal 
coordinates a"" (a = 0, 1), and this induced metric in turn will specify the 
worldsheet measure 

(2) ll7lp/2 

dS = — da° da^ , (99) 

c 

of the timelike 2-surface element spanned by the coordinate variations da° 
and da^ . 

Let us start by supposing that the vorticity distribution is confined 
within a small but finite range specified by displacements 6x~^^ of the 
comoving scalar fields, and then let us take the thin string limit as the size of 
these displacements tends to zero. In this limit it can be seen that the dual 
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vorticity (30) will take the form of a two dimensional Dirac distribution that 



will be expressible in Dirac's notation as an integral over the 2-dimenssional 
world sheet 5^^^ by the formula 

W^""" = -4t7k I W^"'5\x^' - x^'{a]] dS^'^ , (100) 

in which W^^^ is a well behaved antisymmetric tensor on the worldsheet 
representing the concentrated 2 surface vorticity flux. 

In the case of the continuous vorticity distribution considered in the 
preceding section, the vorticity conservation law (^3|) is expressible in its 
dual version in the form 

V^W^""" = . (101) 

When one goes over to the singular limit in which the distribution W^" is 
concentrated in the form (|100|) in the infinitesimal neighbourhood of a single 



flux worldsheet, the conservaton law (101) will translate into a corresponding 



condition on the regular worldsheet supported 2-surface vorticity flux tensor 
W^'^ . Using the abbreviation 

V;. = r]\V, (102) 

for the operator of worldsheet tangential covariant differentiation, where rf^ 
is the "first" fundamental tensor |26, 27 1 of the worldsheet (i.e. the rank-2 



projection operator whose contraction with an abritrary vector at a point 
on the worldsheet projects it onto its tangential part within the worldsheet) 
which will be given by 



r^\ = c~^rP£p,, (103) 

the condition expressing the conservation of the vorticity flux distribution 
in the worldsheet will expressible simply as 

VyW^"' = . (104) 

Since the canonically normalised worldsheet tangential bivector 8^'^ auto- 
matically satifies a conservation condition of the same form 

V^£:''^ = (105) 

as a kinematic identity, it follows that ( |104| ) is interpretable as implying that 
the w 
form 



the worldsheet supported surface vorticity flux tensor W^'^ must have the 



W^"" = kS'"' , (106) 
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where k is a constant on the worldsheet, i.e. V,/K = 0. 



It can be seen that the constant k defined by ( 106 ) will be interpretable 
as the value in the thin string limit of the 2-surface integral of the vortic- 
ity across any small spacelike section through the world tube. By Stoke's 
theorem, using the defining relation (^), this will be equal to the value of 
the Jacobi action around the boundary circuit of the section, which can be 
taken to be any closed curve encircling the vortex string. Thus for any such 
surrounding circuit we shall have 

dS = (f TT^dx" = 2TTnv , (107) 



where v is an integer representing the winding number of the phase 93, i.e. 
v represents the number of individual quantised vortices carrying the flux 
under consideration. This means that if we are considering a string repre- 
senting just a single such quantised vortex, with orientation chosen so that 
= +1 we shall simply have 

K = 2Trn. (108) 

Substituting (^) in formula ( ^5|) it can be seen that one obtains the 
total action integral in the form 

in which the first part is just the 4-dimensional contribution from the irro- 
tational superfiuid outside the vortex, which by (|69|) is 

^(4) = / ^(4) dS''' , C^,, = --^N^'^'N,,, - V{^] , (110) 



while by ( |100| ) and ( |106| ) the second part reduces to a 2-dimensional string 
worldsheet integral given by 

= j dS''' , = ^2^2 _ ^hB^,S>^^ . (Ill) 



The first term in ( 111 ) is like the corresponding action density for a 
simple Nambu Goto (internally structureless) string, as characterised by a 
string worldsheet stress-energy density tensor T'^j^ that is proportional to the 
fundamental tangential projection tensor rj^^^ of the worldsheet as defined by 
( |103| ), except that for an ordinary Nambu-Goto model the effective string 
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tension T is fixed. In the special limit case of the "stiff" fluid model (61) 



which describes the massless axion case, it is well known [|1[ that the vortices 
will have an effective tension T will also be fixed, but for a generic superfluid 
model one obtains 

n = -Tr?",, T=^hH\ (112) 
which shows that the effective tension T will vary proportionally to the 



square of the auxiliary field amplitude, which according to (55) will be given 
on shell by 

$2 = -. (113) 



The second term in ( 111 ) has the form of what is known in the context of 
superstring theory as a "Wess-Zumino" coupling. When such a term turned 
up in the "stiff" massless axion case (|6l|) it was at first interpreted ||l| as a 
"very unusual interaction". However it soon came to be recognised Q as 
a manifestation of an ordinary aerodynamic "lift" type force (arising from 
the Magnus effect) of the kind first evaluated for an aerofoil in the long thin 
(i.e. string type) limit by the Russian theorist Joukowski in the pionneering 
days of subsonic flight a hundred years ago. 



10 The vortex string dynamical equations. 



To work out the effects of the pair of terms in ( |111| ), one must obtain the 
condition for the string type action T^^^ to be invariant with respect to 
an inflnitesimal displacement of the worldsheet generated by an arbitrary 
vector field say. Any such displacement will give rise to corresponding 
"Lagrangian" variations 5$, SB^u, 5g^u (meaning variations as measured 
with respect to local coordinates that are comoving with the displacement) 
of the relevant background fields. Such "Lagrangian" variations will be given 
by corresponding Lie differentiation formulae of the kind introduced above 
for the respective cases of a scalar (0), a closed antisymmetric covariant 
field (p^), and last but in importance not least, the metric itself (19). Thus 
in the trivial case of the scalar $ we shall simply have 

6^ = rv^$ , (114) 

while in the case of the metric we have the familiar expression 

Sg^u = 2V(^/c,) , (115) 
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(which would of course vanish if were not arbitrary but restricted to be a 
Kihing vector generating a background spacetime isometry, i.e. in flatspace 
if A;^ were restricted to be a generator of some Poincare combination of 
translations and Lorentz transformations). Since the antisymmetric field 
Bf^i, is not closed (i.e. since its exterior derivative, the current 3-form N^^p 
is non vavishing) the formula for its Lie derivative is not quite so simple as 
that of its analogue (24) for the vorticity but contains an extra term: the 
relevant formula is 

5B^, = k'^N^.p - 2V[p{B,^pkP) . (116) 

The corresponding variation of the string type action integral specified by 
(|111| ) takes the standard form 

5J(,, = J C^T'^^dgp, + '^W^5Bp, + TS^) ds''^ , (117) 

in which it can be seen, taking account of the variation of the surface element 
( |99| ) due to the induced variation 

'^' = h\h\\v^''^9,u, (118) 



that the relevant surface stress energy momentum tensor T^*^ can be read 
out in the form (|112| ), and that the relevant surface vorticity flux bi- vector 
will evidently be as given by ( |106| ) and ( |108| ), while finally the dynamical 
dual of the amplitude $ can be read out simply as 

y = -T:in^^. (119) 

Since we are only concerned with local variations, there will be no bound- 
ary contribution when we perform the usual operation of integration by parts 
using Green's theorem so as to eliminate gradients of the arbitrary dispace- 
ment vector A;^. The result is thus obtained in the standard form 

5J,,, = j k>^{7p-V,T-^p)dS''' (120) 

in which, after taking account of the surface vorticity fiux conservation law 
( |104| ), the effective surface force density acting in the ensuing dynamical 
equation, 

V,T-^=7^, (121) 



32 



can be read out as 

7^ = |iV/..pW"^+:^V,$. (122) 

In a model of this particular type (as a consequence of the fact that the 
vortex string has no internal degrees of freedom of its own) the tangentially 
projected part of the force balance equation ( |12lD , i.e. the part that is ob- 
tained by contracting it with the fundamental tensor rj^^p given by ( [L03D , will 
automatically be satisfied as a kinematic identity. The only dynamical infor- 
mation in (121) is the part obtained by contracting it with the complement 



of the tangential projection tensor r/'^p, namely the orthogonal projection 
tensor defined by 

1-^^=9%- tp- (123) 

It can be seen that the remaining, purely orthogonal part of the left 
hand side of ( |121| ) can be evaluated in terms of the extrinsic geometry of 
the imbedding, without knowledge of the gradient of T^^^ , using the fact that, 
since the surface stress momentum energy tensor T^^ is purely tangential, 
i.e. since its contraction with J.'^ vanishes, an integration by parts leads to 
the identity 

IF^ V.T'^M = T^-i^^/ , (124) 

where the K^yP is the second fundamental tensor as defined [26| in terms of 
tangential differentiation of the first fundamental tensor rjP^ by 

K^,P = r^^,V^r^>'„. (125) 



One thus obtains an expression of the standard form |27] 

T^^if^/ =±^^ 7^ , (126) 
for the extrinsic part of the dynamical equation (121), i.e. the part that 



governs the evolution of the worldsheet itself, which is the only part there 
is in the present case. 

Due to the Nambu-Goto like form (|112| ) obtained for T^^ in the simple 
vortex model under consideration here, it will not be necessary to work out 
the complete second fundamental tensor K^^'^ in the present case, but only 
its trace, the curvature vector 

RP = K\P = Vyi'P , (127) 



in terms of which it can be seen that the equation of motion ( [12 61) will reduce 
to the form 

TKp =Ay v,T + TTne^'N^yp , (128) 
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with the effective tension T given by ( |112| ), according to which it is pro- 
portional to the square of the auxihary field <I> in the ambient superfluid 
background. The first term on the right of ( |128| ), allowing for the effect of 
any non-uniformity of this field, was not needed in the special limit of the 
"stiff" fluid model ( |6l] ) characterised by a fixed value of which describes 
the massless axion case to which earlier studies of relativistic vortex string 
dynamics |||, Q were restricted. 

11 Vorton equilibrium in stationary background 

As a simple exercise, let us apply the foregoing theory to the problem of a 
vorton, i.e. an equilibrium state of a vortex ring in a uniform background. 
As in the well known non relativistic version of this problem, relative motion 
is needed to provide the Joukowsky type force (due to the Magnus effect) 
that supports the string against its own tension. A superfluid vorton differs 
from local cosmic vorton in that the latter is supported not by a 

Magnus effect, but by the centrifugal effect of circulating current. 

By definition, an equilibrium state is invariant with respect to the action 
of a time displacement vector, say, which, to preserve the background 
metric must satisfy the Killing equation given according to (|l9|) by 2V(^A:jy) = 
0. This holds as a triviality for an ordinary flat space time displacement 
vector, which satisfies 

Vf,K = . (129) 

Invariance requires that the vortex worldsheet be tangent to this vector, 
whose components in a corresponding Minkowski coordinate system will be 
given by /c^ <-> {1, 0, 0, 0} with normalisation k^^k^ = —(? . 

The condition that the superfluid background is not just stationary but 
uniform with velocity v relative to such a frame means that for suitably 
aligned space axes its 4-velocity will be given by <-> 7{1, 0, 0} with 
7 = (1 — u^/c^)"^/^. The assumption that the ambient fluid is uniform 
means that the gradient term will drop out of the dynamical equation ( |128D , 
with the consequence that the background number density n will also cancel 
out, leaving an equilibrium condition that is expressible just in terms of the 
4-momentum covector tx^ = fiUy in the form 

niKp = 2S^'''e^,p„^:'' . (130) 

In such a stationary case it is easy to evaluate the tangent bivector 
8^^ and the worldsheet curvature vector in terms of the unit tangent 
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vector, say, that is orthogonal to within the worldsheet (which to be 
uniquely specified requires a choice of orientation). In terms of this vector, 
as characterised by e^/c^ = 0, e'^e^ = 1, we shall simply have 

and (independently of the choice of orientation) the corresponding funda- 
mental tensor will be given by 

r/^'^ = -^M^^i'^ + e^e^ (132) 



It can be seen from the defining relations ( 102 ) and ( |127D that when ( |129| ) 
is satisfied we shall simply have 

K^' = e'^V^e'^ . (133) 

The vorton equilibrium states in which we are interested will be symmet- 
ric about an axis aligned with the relative flow, which we take to be the 
axis of our Minkowski coordinate system. The worldsheet will be circular, 
with fixed radius r say, meaning that any point on it will be located relative 
to the axis by a radial displacment vector that can be parametrised in terms 
of the corresponding axial angle 9 in the form r{0, 0, ,cos0, sin0} 

with r'^r^ = r^. The corresponding unit tangent vector will be given by 
e'^ <-> {0,0, , —sin 0, cos 0} and thus the corresponding curvature vector is 
easily seen to be given by = r^'^r^ r^^{0, 0, ,cos0, sin0}. 

From this familiar result, namely that the curvature vector K'^ of a circle 
of radius r is inwardly directed with magnitude r~^, it can be seen that the 
solution to the equilibrium condition ( |130| ) will be obtained for a vorton 
radius given by 

%l , ^ 

in which it is to be recalled that, according to ([T^), I is an order of unity 
factor that can be estimated as the (natural) logarithm of the ratio of the 
vorton area to the sectional area of the vortex core. 

The vagueness of the prescription for the / symptomises an inherent lim- 
itation of any attempt (whether in a relativistic or in a Newtonian frame- 
work) to treat a "global" vortex as if it were a locally confined string type 
phenomenon, despite the "infra-red divergence" of its energy and tension, 
which can only be made finite by a long range cut off. This problem does 
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not arise for the vortex defects in an ordinary metallic type II supercon- 
ductor, in which the local gauge coupling to the electromagnetic field en- 
sures an exponential fall off ensuring convergence without any need for a 
cut off. In the cosmic string context, there are locally gauged examples for 
which a string type description is highly accurate there are also examples 
of non-local vortex defects, such as the axionic case ^ (and cases involving 
electromagnetic [^] or gravitational coupling) for which the thin string 
limit description is less satisfactory. If the vortex core radius 6 is relatively 
small, the problem of sensitivity to the long range cut off A is mitigated by 
the fact that the dependence is only logarithmic. Thus a formula such as 
( |134| ) (or its well known non-relativistic anologue, in which the 7 factor is 
omitted, and /j, is replaced by a fixed mass m) can be usefully applied to 
macroscopic vorton configurations, but becomes quantitatively unreliable for 
describing the microscopic vortons, known as "rotons", that are important 
in the analysis of thermal excitations. 

Whereas the problem of specification of the cut off can limit the util- 
ity of the string type description (111) for application to the dynamics of 
individual vortices, this reservation does not apply to the macroscopically 
averaged description for an Abrikosov type lattice of aligned vortices as given 
by (^), which is on a much sounder footing: in this case the relevant cut off 
is unambiguously determined by the lattice spacing. Although its accuracy 
is more questionable, the single string picture provides useful insight into 
the intepretation of the - at first sight rather mysterious - dynamical equa- 
tion (|82|) for the more robust macroscopically averaged model, which is in 
fact formally identical to the more intuitively interpretable string dynamical 
equations ( |128D . 
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